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Abstract: 

For the Ising model, the spin magnetization transition is equivalent to the percolation 
transition of Fortuin-Kasteleyn clusters; this result remains valid also for the conven- 
tional continuous spin Ising model. The investigation of more general continuous spin 
models may help to obtain a percolation formulation for the critical behaviour in SU{2) 
gauge theory. We therefore study a broad class of theories, introducing spin distribution 
functions, longer range interactions and self-interaction terms. The thermal behaviour of 
each model turns out to be in the Ising universality class. The corresponding percolation 
formulations are then obtained by extending the Fortuin-Kasteleyn cluster definition; in 
several cases they illustrate recent rigorous results. 



Introduction 

Percolation theory 0] has been successfully applied to describe continuous thermal 
phase transitions as geometrical transitions: at the critical threshold, suitably defined 
clusters become infinite (see 0-0). It is moreover possible to establish a correspondence 
between thermal and cluster variables, such that the critical exponents of corresponding 
variables coincide Q]. This provides an intuitive geometric description of the mechanism 
of the thermal transition; e.g., in the case of the Ising model, we can identify clusters 
as magnetic domains, so that magnetization sets in when they span the lattice. In a 
recent study 0, we have tried to extend this picture to the confinement-deconfinement 
transition of SU{2) gauge theory. This theory leads to a continuous phase transition 
which belongs to the universality class of the Ising model, due to the common Z{2) global 
symmetry of the corresponding Hamiltonians However, our result is valid only in the 
strong coupling limit of SU{2), where the partition function can be well approximated 
by that of a continuous spin Ising model. In this model, the correspondence between 
percolation and thermal variables can be rigorously established |jTD|, rT|. For a general 



study of SU{2) gauge theory, the weak coupling limit must be taken into account, and 
there the approximations used in |p are no longer valid. 

To address the deconfinement problem more generally, it may be helpful to look for 
an effective theory for SU{2). An interesting attempt in this direction was proposed 
some time ago in the search for the fixed point of the theory by means of block-spin 
transformations [Q. There the original Polyakov loop configurations were projected onto 



configurations of Ising-like spins according to the sign of the Polyakov loop at each lattice 



site. This assumes that for the criticahty only the global Z{2) symmetry of the theory is 
important. As a result of this projection one obtains an effective Hamiltonian with only 
short-range couplings, although some of the resulting couplings correspond to interactions 
beyond the fundamental nearest-neighbour form. Moreover, the assumption that only the 
Z{2) symmetry of the theory is relevant must be tested. 

In this work we investigate several spin models in order to get a better understanding 
of the behaviour of the SU{2) effective theory. Our investigations are based on lattice 
Monte Carlo simulations of the various models. 

We focus on continuous spin Ising models, in which the individual spins Sj at each 
lattice site can take on all values in the finite range [—1,1], with the distribution of the 
spins governed by a spin distribution function /(sj). For /(sj) = — 1) V i and 

only nearest-neighbour (NN) interactions, we recover the conventional Ising model. For 
f{si) = 1 \f i and only NN interactions, we obtain the continuous spin model studied in 
[p!0| . Here we will consider four more general models of this type {d denotes the space 
dimension) : 

A) d = 2, only NN interactions, but a non-uniform spin distribution f{si); 

B) d = 2, f{si) = 1 V i, NN and diagonal next-nearest-neighbour (NTN) interaction (Fig. 

1); 

C) d = 3, f{si) = 1 V i, NN and two types of NTN interactions (see Fig. 2); 

D) as case C), but including an additional self- interaction term proportional to V i. 
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Figure 1. Spin interactions in 
model B. 








\ ^ 

\ 


. J3IAG. 


\ 

V 




NTN ^ 












If' 





Figure 2. Spin interactions in 
model C. 



We will first present a detailed numerical investigation of the thermal critical behaviour 
and show that each model belongs to the same universality class as the Ising model. 
Subsequently, we provide a suitable cluster definition in order to interpret the thermal 
phase transition of the models as a geometrical percolation transition. Concerning the 
second issue, it was recently proven rigorously that the result of [|1^ can be extended 
to the models A, B and C |TT|. We believe, however, that it is useful to illustrate it by 



exploiting the power of the finite-size scaling analysis. On the other hand, model D has 
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not been investigated so far; as self-interactions play an important role in gauge theories, 
it is interesting to examine what happens in this case. 



Model A: Spin Distribution Function 

The spin variable s can take all values in [—1,1]. It is convenient to rewrite it in the 
form 

s = sign{s)a, (1) 

separating the sign from the amplitude of the spin. Since we want to retain the Z(2) 
symmetry of the Hamiltonian, the signs of the spins are equidistributed. In contrast, 
the amplitudes can in general be weighted in different ways by choosing a distribution 
function /(cr). The partition function of such a model then has the form 

^(^)=n/ daif{a^) exp{Kj2siSj}, (2) 

where k = J/T, with J denoting the spin coupling energy and T the temperature. Since 
we want to study models with continuous transitions, the choice of the function /(a) is 
not arbitrary. It can be proved that it must obey certain regularity conditions, which 
are not very restrictive, however |jl3|. In [|l^, the amplitudes were equidistributed, with 



/((t) = 1. Here we consider the following form for /(cr): 



/(a) = Vl^ (3) 

which is the Haar measure of the SU(2) group; it appears in the strong coupling ap- 
proximation of SU(2) 0. To define a cluster, we use exactly the same procedure as in 
p!0| |, i.e., we join nearest-neighbouring spins of the same sign with the bond probability 



p = 1 — exp{—2Kaiaj), where cxj and aj are the amplitudes of the two spins. We use 
free boundary conditions and define percolation when a cluster spans the lattice in each 
of the two directions. This is done to exclude the possibility that, due to the finite size 
of our lattices, one could find more than one percolation cluster, which would make the 
definition of our variables ambiguous. We shall use this definition for all the models to 
be considered hereQ. 

The percolation variables to be studied are: the percolation strength P, 

^ size of the percolating cluster 
volume of the lattice 

and the average cluster size S, 

S = %^. (5) 

Here rig is the number of clusters of size s and the sums exclude the percolating cluster. 
Besides these two, we evaluate another very useful variable. For a given lattice size and 
a value of k, we determine the relative fraction of configurations containing a percolating 

^In three dimensions even this definition of a spanning cluster does not exclude the possibility of 
having more than one such cluster in the same configuration. Nevertheless the occurrence of these cases 
is so small that we can safely ignore them. 
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cluster. This variable, that we shall call percolation cumulant and indicate as 7r(K), is a 
scaling function, analogous to the Binder cumulant gr{i^) [ H - 

In addition to these percolation variables, we also measure the lattice averages of 
energy and magnetization, in order to study the thermal transition as well. From the 
magnetization M we calculate in particular the physical susceptibility Xy 

X{n) = V ((M^) - (Mf) (6) 

and the Binder cumulant gr, 



9r{^) = ^ ■ (7) 



We have performed simulations on four lattice sizes, from 64^ to 240^. Our algorithm 
consists in heatbath steps for the update of the spin amplitudes followed by Wolff-like 
cluster updates for the flipping of the signs. That is basically the same method as used 
in ||10|, although the heatbath procedure is slightly modified to take into account the 



presence of the distribution function /(a). The proof that the detailed balance condition 



is satisfied for such a cluster update is obtained following the derivations in [0] - [jTH 
Flipping the spins reduces drastically the autocorrelation time, which allowed us to collect 
extensive data and thus reduce the errors. 

Fig. 3 shows a comparison between the thermal variable gr{i^) and the percolation 
cumulant 7r(/t), both as functions of the temperature variable k, for different lattice sizes. 
The points where the lines cross are the infinite volume thresholds for the thermal and the 
geometrical transition, respectively. The agreement between the two points is remarkable. 

Since the percolation cumulant is a scaling function, we can already get an indication 
about the critical exponents of the percolation transition. Given the critical point and the 
exponent z/, a rescaling of the percolation cumulant as a function of {k — K(.r)L^^'^ should 
give us the same function for each lattice size (L is the linear dimension of the lattice). 
Figs. 4 and 5 show the rescaled percolation cumulant, using Hcr = 1-3888 with the random 
percolation exponent z/ = 4/3 and Ising exponent z^ = 1, respectively. The figures clearly 
show scaling for the Ising exponent and no scaling for the random percolation exponent. 

The exponent f3 governs percolation strength and magnetization, while 7 governs 
cluster size and susceptibility. To determine them, we have performed high-statistics 
simulations around the critical point, with the number of measurements for each value 
of the coupling varying from 50000 to 100000. We have then used the method |9| to 



determine the values of the exponents. The results are shown in Table I, where the first 
row shows the calculated percolation results, the second the corresponding magnetization 
calculations, and the third the analytically known Ising model values. It is evident that the 
percolation behaviour coincides fully with the thermal critical behaviour. This conclusion 



holds in general for all admissable distribution functions, as it was shown in |TT 

Model B: Next-to-Nearest Neighbour Interactions 

We now want to study how to correctly define clusters in the presence of more than 
the standard nearest-neighbour interactions. We have now two terms, with a Hamiltonian 
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of the form 




(8) 



{hi) {hi) 



where the first sum describes nearest-neighbour and the second diagonal next-to-nearest 
neighbour interactions (Fig. 1). Since longer range interactions are generally weaker, we 
fix the ratio between the two couplings at Jnn/Jntn = 10; however, we do not believe 
that our results depend on the choice of the couplings, as long as both are ferromagnetic. 
To define clusters, we now extend the Conigho-Klein method and define for each two spins 
i,j of the same sign, for NN as well as NTN, a bond probability 



where x specifies ktvat = Jnn/T and kntn = Jntn/T, respectively. We recall that in 
model B, as well as in C and D, we assume a spin equidistribution, i.e., /(ex) = 1. 

We have studied model B using two different Monte Carlo algorithms, in order to test 
if a Woffi-typc algorithm can also be applied in the presence of NTN interactions. The 
first is the standard Metropolis update, while the second alternates heat bath steps and 
a generahzed Wolff hipping, for which the clusters are formed taking into account both 
interactions. The generalization of the cluster update is trivial. After several runs, some 
with high statistics, we found excellent agreement with the Metropolis results in all cases. 
So, the mixed algorithm with heat bath and Wolff flippings appears to remain viable also 
in the presence of more than the standard NN interaction. Subsequently we have therefore 
used this mixed algorithm. The lattice sizes for model B ranged from 100^ to 400^. 

As for model A, we now show for model B the comparison between percolation cu- 
mulant and Binder cumulant, in order to test that the critical points coincide (Fig. 6). 
We then again rescale the percolation cumulant, using the critical k determined in Fig. 
6 together with the exponent u from the d = 2 random percolation transition (z/ = 4/3) 
and the d = 2 Ising model (t/ = 1). Figs. 7 and 8 indicate that again the correct exponent 
is that of the Ising model. 

The determination of the other two exponents /? and 7 by means of the method 
confirms that indeed both the thermal and the geometrical exponents belong to the Ising 
universality class (Table II). 

Model C: Extension to Three Dimensions 

We have here repeated the study for a d = 3 model with three different interactions 
(Fig. 2). To fix the model, we have to specify the ratios of the nearest-neighbour coupling 
Katat to Untn and Kdiag- We chose them to be 10:2 and 10:1, respectively. Our d — 3 
calculations are performed on lattices ranging from 12^ to 40^. 

Also here we have first compared the results from a mixed algorithm of the same 
kind as for the previous case to tliosc from a standard MetropoUs algorithm; again, the 
agreement is very good. Figs. 9, 10 and 11 then show the comparison of the thresholds 
and the scaling of the percolation cumulant. As before, the correspondence between 
percolation and thermal variables is evident (Table III). 



= 1 - exp{-2n^aiaj}, 



(9) 
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Model D: Adding Self-Interaction 

From what we have seen up to now, it seems to be clear that the correct cluster 
definition can readily be extended to models including several (ferromagnetic) spin-spin 
interactions. However, such terms are not the only possible interactions in a model with 
Z(2) symmetry and a continuous transition. There could be anti-ferromagnetic spin-spin 
couplings as well as multispin terms, coupling an even number of spins greater than two 
(four, six, etc.). Moreover, since the spins are continuous, the presence of self-interaction 
terms is possible, determined by s^, s^, etc. The treatment for antiferromagnetic and 
multispin couplings so far remains an open question. In contrast, self-interactions are not 
expected to play a role in the cluster building, since such terms do not relate different 
spins. Therefore, we test a cluster definition ignoring any self-interaction term. 

We thus consider in Model D the same interactions as in Model C, plus a term propor- 
tional to JoSi-sf. We chose a negative value for the self-interaction coupling Jq; this is 
the more interesting case since the corresponding interaction tries to resist the approach 
of the system to the ground state at low temperatures {a = 1 everywhere). The ratios of 
the NN coupling to the others were chosen as Jatat : Jntn '■ Jdiag :|Jo|=6:2:l:2. 

Again we first verify the viability of the mixed algorithm and then determine the 
rescaled percolation cumulant (Fig. 12) and the critical exponents (Table IV). It is evident 
that the percolation and the thermal transition again fall into the Ising universality class. 

Conclusions 



We have performed a complete investigation of several continuous spin models, ob- 
tained from the simple model studied in [10| by introducing a distribution function for 



the spin amplitudes, longer-ranged interactions and self-interactions. The thermal be- 
haviour leads in all cases to the Ising universality class. All systems admit equivalent 
percolation formulations; for models A, B and C, we illustrate numerically the analytical 
results found in [|TT|. All the features of the models we studied are expected to enter in 



the further extension to SU(2) gauge theories, and the fact that they do not interfere with 
a percolation formulation of critical behaviour provides support for such efforts. 
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Table I. Thermal and percolation thresholds and exponents of Model A. 
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Table II. Thermal and percolation thresholds and exponents of Model B. 
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Table III. Thermal and percolation thresholds and exponents of Model C. 
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Table IV. Thermal and percolation thresholds and exponents of Model D. 
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Figure 3. Comparison of the thermal and the geometrical critical point for 
Model A obtained respectively from the Binder cumulant gr and the percolation 
cumulant. 
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Figure 4. Rescalod percolation cumulants for Model A using the 2-dimensional 
random percolation exponent v = 4/3. 
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Figure 5. Rescaled percolation cumulants for Model A using the 2-dimensional 
Ising exponent v = 1. 
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Figure 6. Comparison of the thermal and the geometrical critical point for 
Model B obtained respectively from the Binder cumulant gr and the percolation 
cumulant. 
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Figure 7. Rcscalcd percolation cumulants for Model B using the 2-dimensional 
random percolation exponent v = 4/3. 
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Figure 8. Rescaled percolation cumulants for Model B using the 2-dimensional 
Ising exponent v = 1. 
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Figure 9. Comparison of the thermal and the geometrical critical point for 
Model C obtained respectively from the Binder cumulant gr and the percolation 
cumulant. 
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Figure 10. Rcscalod percolation cumulants for Model C using the 3-dimensional 
random percolation exponent = 0.88. 



E 

o 



0.8 
0.6 
0.4 
0.2 



size 12, 

size 20, X 

size 30, ^ - X - 

size40'^ > Q 



X X- 







-2.5 



-1.5 



-1 



-0.5 

I (l/v) 



0.5 1 



1.5 



Figure 11. Rescaled percolation cumulants for Model C using the 3-dimensional 
Ising exponent v = 0.6289. 
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Figure 12. Rescaled percolation cumulants for Model D using the 3-diniensional 
Ising exponent v = 0.6289. 
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